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ABSTRACT 

Using the operator product expansion and nonrelativistic QCD we study the 
production of charmonium hybrids in B decays. We express the decay rate in 
terms of a few matrix elements which eventually will be fixed by experimental 
measurements or calculated on the lattice. While the Fock state expansion is 
problematic for hybrids, the operator product expansion still provides a model 
independent framework for phenomenological calculations of hybrid production 
and decay. We then use a simple flux tube model to estimate the branching ratio 
B ij)g + X, where ik g is a J PC = Cd hybrid, the large production of which 
could help resolve the low charm multiplicity observed in B decays. 


The hybrid states remain an enigmatic feature of hadronic physics. Their existence 
is predicted in QCD, yet most of their properties have been studied only within models 
or on the lattice. This is true even of quarkonium hybrids, because, unlike ordinary 
quarkonium, hybrids are characterized by excitations of the nonperturbative gluon field. 
Even for large mass rriQ, there is no reason for a QQ hybrid to be a compact object, and 
neither the lattice nor models support such a picture. Contrary to the hybrids, in the 
limit mQ —> oo ordinary quarkonium admits a controlled Fock state decomposition, W in 
which the leading term is a QQ in a color singlet with fixed quantum numbers 2S+1 Lj. 
Higher order terms, such as color octet configurations, are suppressed by powers of 
a s (m Q ) v, where v is the relative velocity of the Q and the Q. The expansion is 
possible because for large m<q the quarkonium is a small state, of size (m^w) -1 , whose 
interaction with the color field is governed by a multipole expansion. This velocity 
power counting can be made explicit by a suitable rescaling of the fields in the effective 
Lagrangian of NRQCD .□ By contrast, the minimum of the hybrid meson potential is at 
a separation r 0 of order 1/Aq CD , independent of rriq. For large m,Q, fluctuations about 
r 0 are small, but the state itself is not compact. This is a generic feature of hybrid 
models; for example, in a constituent gluon model, the Q and Q are in a color octet 
and repel each other at short distances. At large tuq, one expects a situation somewhat 
like a heavy rigid rotor, with nearly degenerate rotational bands of states. This behavior 
has been observed in lattice studies.B We can make the argument explicit by rescaling 
the coordinates in the NRQCD Lagrangian by their “natural” size, x = A^X, t = A t T, 


where, contrary to the quarkonium case! where A^ = l/(mv) we have X x = 1/A qcd- 
This rescaling explicitly takes into account the fact that in the limit mg —> oo hybrid 
quarkonium state remains finite. A NRQCD Lagrangian can be written as 
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It implies that <9° is of the same order as V 2 /2 m Q , leading to X t = rriQ\ 2 x = mq/A q C d . 
The helds in ([[]) should also be rescaled as Q — AqT, A 0 = A^oA 0 , A = A^A 
where from the fact that the change of integration variables brings another factor 
of A 3 A 4 we obtain that A q = A/. The behavior of the gauge part of the Lagrangian 
is slightly more subtle, but it is easy to convince yourself that in the heavy 
quark limit m Q /A QC d —> oo the terms (V x A) 2 —> 7tiqA^Aa(V x A) 2 and (VA 0 ) 2 —> 
m Q XlX A o(yA 0 ) 2 are the leading ones with Aa = A a o = \j {jnqXl)- 1 = A QCD \/A QC D/m Q 
This gives the rescaled NRQCD Lagrangian for hybrid quarkonia 
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Clearly, all the helds and derivatives in f^j are dimensionless, so all suppression factors 
explicitly appear in front of the relevant operators providing a power counting scheme of 
the theory. Note that Eq.(^j) implies that the Coulomb (i.e. A 0 ) gluon interactions must 
be resummed to all orders. Hereafter we shall use conventional helds and derivatives 
keeping Eq.(pj) in mind to maintain correct power counting of the relevant operators. 

The absence of a controlled Fock state expansion for hybrid quarkonium can be 
understood as follows. For the ordinary heavy quarkonium, a soft glue configuration 
has a small overlap with the compact two-quark state, since the Compton wavelength 
of the soft gluons, of the order of A qq Dj is much larger than the distance between the 
quarks. By contrast, the size of the heavy hybrid state is finite, of order A^ CD in the 
heavy quark limit, allowing for a significant interaction with the soft gluonic modes. 

Nevertheless, there is still a sort of Fock state expansion for quarkonium hybrids. 
Although the gluonic degrees of freedom are intrinsically nonperturbative, we can label 
the various components according to the quantum numbers of the QQ. For example, 
for the 0 +_ charmonium hybrid state, we can write a decomposition of the form 

IVb) = A|cc( 3 S'i)i , 8 + g±) + H|cc( 3 P/)i , 8 + # 2 ) + (- |cc( 1 S'o)i,8 + 9z) + • • • , (3) 


where the subscripts 1, 8 indicate the color of the cc pair and the g* are various gluonic 
configurations. While one is always free to do this, here there is no model-independent 
hierarchy among the coefficients. This is related to the lack of a unique cc “baseline” 
state, and the lack of a multipole expansion to govern transitions between the different 
Fock components. Nonetheless, within particular models it is typically the case that a 
given hybrid is dominated by one or two cc configurations. Hence it is useful to organize 
our calculation of hybrid production in B decays with a decomposition such as (^j) . 















As an application of the formalism of nonrelativistic QCD to the physics of 
quarkonium hybrids, we calculate hybrid charmonium production in B decayshl As 
usual, we first calculate the production rate in the framework of full QCD, separating 
hard and soft degrees of freedom. Then we compute the rate in NRQCD, and use the 
perturbative calculation to match the coefficients of the NRQCD matrix elements. The 
inclusive B meson decay rate to a hybrid charmonium state B —> + X is 
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where P Sc = \cc + X 3 ) (cc + X 3 | is a projection operator onto the cc state in the ip g , 
P = sT^T a b. The Eq. (£|) is then matched with the equivalent expression in NRQCD, 
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where O^ 9 represents a set of NRQCD operators of increasing dimension d n 
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Here P^ g is a projection onto the hybrid state | tp g ) which is normalized nonrelativisti- 
cally. Note that each derivative insertion brings a power of the velocity of the heavy 
quark, p/m c ~ v: Og 9 ( 3 Pi) is suppressed compared to Og 9 ( 3 S'i) by v 2 . We are unable 
to calculate matrix elements of these operators in a model independent fashion, so as 
usual they are left as free parameters. Since the same matrix elements would govern 
the total annihilation width of the hybrid or the hybrid photoproduction cross section, 
one should extract the leading matrix elements from a small number of experiments. 

Performing the perturbative calculation of , boosting the relativistic four- 
component spinors u c to the cc center of mass frame and then replacing them by 
nonrelativistic two-component spinors £ and rj we arrive at a simple expression 
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where /7 C = m 2 /ml and ^ = rrvfy/m 2 . Eq. (0) expresses the decay rate in terms of 
matrix elements of operators of increasing dimension. 

To perform the matching, we have to calculate T(B —> ifjg + A") in perturbative 
NRQCD using Eqs.(j3j) and (p|) and then evaluate the NRQCD coefficients by matching 
full NRQCD states to perturbative ones. The partial decay rate for B —> xjj g + X is 
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where we have deffiied T 0 = Cl\V c bV* s \ 2 G 2 F m\/lAA'K. Since the Fock space decomposi¬ 
tion of hybrid states is not controlled by a multipole expansion, the number of matrix 
elements important for production processes can be reduced by keeping those that are 
leading in powers of l/m c . The formula (||) determines the decay rate of B into a given 
hybrid charmonium state in terms of universal parameters which must be fixed by 
other experiments. Using the flux tube model! to estimate the matrix elements ([]) we 
End! 

T(B i/j g + X) < (1 x 10 -3 ) T(b —7 ccs). (9) 


Thus we estimate that the production of this hybrid charmonium is more than a factor 
of ten lower than that of typical ordinary charmonium states. It is expected that the 
produced hybrids, although heavy, would not primarily decay into the final states 
containing charmed quarksj2l thus offering a solution to the problem of low charm 
multiplicity in B decays. However, given the tiny overall magnitude of the result, we 
conclude that it is quite unlikely that this mechanism could play an important role in 
“hiding” the charm produced in b —*■ ccs transitions. It however places the studies of 
hybrid charmonia well within the reach of the incoming B factories.Q 
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